##########################################
####### MAPPING THE CZECH ELECTIONS ######
##########################################

# DATA AND ANALYSIS  
# Jakub Lysek, Tomáš Lebeda, Jiří Pánek


# import datasetu
setwd("~/ArcCR500/Výsledky voleb/Sněmovní volby - 1996 - 2017")

library(haven)
PS_KOMPLET_17_96 <- read_sav("~/ArcCR500/Výsledky voleb/Sněmovní volby - 1996 - 2017/PS_KOMPLET_17_96.sav")
#View(PS_KOMPLET_17_96) #POZOR data jsou bez Prahy! pro analýzy sněmovních voleb 2017 a 2013 je vhodnější použít druhý soubor
FINAL_17 <- read_sav("~/ArcCR500/Výsledky voleb/Sněmovní volby - 1996 - 2017/FINAL_17.sav")



#Standardizace vybraných proměnných - dataset PS_KOMPLET_17_96
PS_KOMPLET_17_96$sd_PAR_2017_4_Českástr.sociálnědemokrat._p <- scale(PS_KOMPLET_17_96$PAR_2017_4_Českástr.sociálnědemokrat._p)
PS_KOMPLET_17_96$sd_PAR13_1_Česká_strsociálně_demokrat_p <- scale(PS_KOMPLET_17_96$PAR13_1_Česká_strsociálně_demokrat_p)
PS_KOMPLET_17_96$sd_PAR10_9_Česká_strsociálně_demokrat_p <- scale(PS_KOMPLET_17_96$PAR10_9_Česká_strsociálně_demokrat_p)
PS_KOMPLET_17_96$sd_PAR06_10_Česká_strsociálně_demokrat_p <- scale(PS_KOMPLET_17_96$PAR06_10_Česká_strsociálně_demokrat_p)
PS_KOMPLET_17_96$sd_PAR02_3_Česká_strsociálně_demokrat_p <- scale(PS_KOMPLET_17_96$PAR02_3_Česká_strsociálně_demokrat_p)
PS_KOMPLET_17_96$sd_PAR98_7_Česká_strsociálně_demokrat_p <- scale(PS_KOMPLET_17_96$PAR98_7_Česká_strsociálně_demokrat_p)
PS_KOMPLET_17_96$sd_PAR96_2_Česká_strsociáldemokratická_p <- scale(PS_KOMPLET_17_96$PAR96_2_Česká_strsociáldemokratická_p)

PS_KOMPLET_17_96$sd_PAR_2017_8_Komunistickástr.ČechaMoravy_p <- scale(PS_KOMPLET_17_96$PAR_2017_8_Komunistickástr.ČechaMoravy_p)
PS_KOMPLET_17_96$sd_PAR13_21_Komunistická_strČech_a_Moravy_p <- scale(PS_KOMPLET_17_96$PAR13_21_Komunistická_strČech_a_Moravy_p)
PS_KOMPLET_17_96$sd_PAR10_6_Komunistická_strČech_a_Moravy_p <- scale(PS_KOMPLET_17_96$PAR10_6_Komunistická_strČech_a_Moravy_p)
PS_KOMPLET_17_96$sd_PAR06_20_Komunistická_strČech_a_Moravy_p <- scale(PS_KOMPLET_17_96$PAR06_20_Komunistická_strČech_a_Moravy_p)
PS_KOMPLET_17_96$sd_PAR02_23_Komunistická_strČech_a_Moravy_p <- scale(PS_KOMPLET_17_96$PAR02_23_Komunistická_strČech_a_Moravy_p)
PS_KOMPLET_17_96$sd_PAR98_9_Komunistická_strČech_a_Moravy_p <- scale(PS_KOMPLET_17_96$PAR98_9_Komunistická_strČech_a_Moravy_p)
PS_KOMPLET_17_96$sd_PAR96_15_Komunistická_strČech_a_Moravy_p <- scale(PS_KOMPLET_17_96$PAR96_15_Komunistická_strČech_a_Moravy_p)

#Výpočet nových proměnných
PS_KOMPLET_17_96$ROZ_CSSD_98_96 <- (PS_KOMPLET_17_96$PAR98_7_Česká_strsociálně_demokrat_p - PS_KOMPLET_17_96$PAR96_2_Česká_strsociáldemokratická_p)
PS_KOMPLET_17_96$ROZ_CSSD_02_98 <- (PS_KOMPLET_17_96$PAR02_3_Česká_strsociálně_demokrat_p - PS_KOMPLET_17_96$PAR98_7_Česká_strsociálně_demokrat_p)
PS_KOMPLET_17_96$ROZ_CSSD_06_02 <- (PS_KOMPLET_17_96$PAR06_10_Česká_strsociálně_demokrat_p - PS_KOMPLET_17_96$PAR02_3_Česká_strsociálně_demokrat_p)
PS_KOMPLET_17_96$ROZ_CSSD_10_06 <- (PS_KOMPLET_17_96$PAR10_9_Česká_strsociálně_demokrat_p - PS_KOMPLET_17_96$PAR06_10_Česká_strsociálně_demokrat_p)
PS_KOMPLET_17_96$ROZ_CSSD_13_10 <- (PS_KOMPLET_17_96$PAR13_1_Česká_strsociálně_demokrat_p - PS_KOMPLET_17_96$PAR10_9_Česká_strsociálně_demokrat_p)
PS_KOMPLET_17_96$ROZ_CSSD_17_13 <- (PS_KOMPLET_17_96$PAR_2017_4_Českástr.sociálnědemokrat._p - PS_KOMPLET_17_96$PAR13_1_Česká_strsociálně_demokrat_p)

PS_KOMPLET_17_96$ROZ_KSCM_98_96 <- (PS_KOMPLET_17_96$PAR98_9_Komunistická_strČech_a_Moravy_p - PS_KOMPLET_17_96$PAR96_15_Komunistická_strČech_a_Moravy_p)
PS_KOMPLET_17_96$ROZ_KSCM_02_98 <- (PS_KOMPLET_17_96$PAR02_23_Komunistická_strČech_a_Moravy_p - PS_KOMPLET_17_96$PAR98_9_Komunistická_strČech_a_Moravy_p)
PS_KOMPLET_17_96$ROZ_KSCM_06_02 <- (PS_KOMPLET_17_96$PAR06_20_Komunistická_strČech_a_Moravy_p - PS_KOMPLET_17_96$PAR02_23_Komunistická_strČech_a_Moravy_p)
PS_KOMPLET_17_96$ROZ_KSCM_10_06 <- (PS_KOMPLET_17_96$PAR10_6_Komunistická_strČech_a_Moravy_p - PS_KOMPLET_17_96$PAR06_20_Komunistická_strČech_a_Moravy_p)
PS_KOMPLET_17_96$ROZ_KSCM_13_10 <- (PS_KOMPLET_17_96$PAR13_21_Komunistická_strČech_a_Moravy_p - PS_KOMPLET_17_96$PAR10_6_Komunistická_strČech_a_Moravy_p)
PS_KOMPLET_17_96$ROZ_KSCM_17_13 <- (PS_KOMPLET_17_96$PAR_2017_8_Komunistickástr.ČechaMoravy_p - PS_KOMPLET_17_96$PAR13_21_Komunistická_strČech_a_Moravy_p)

PS_KOMPLET_17_96$sd_INHABITANTS_2018_LOGGED <- scale(PS_KOMPLET_17_96$INHABITANTS_2018_LOGGED)
PS_KOMPLET_17_96$sd_EXEKUCE_OBCE <- scale(PS_KOMPLET_17_96$EXEKUCE_OBCE)


#zaokrouhlení na dvě místa
PS_KOMPLET_17_96[,'PAR13_1_Česká_strsociálně_demokrat_p']=round(PS_KOMPLET_17_96[,'PAR13_1_Česká_strsociálně_demokrat_p'],2)
PS_KOMPLET_17_96[,'PAR_2017_4_Českástr.sociálnědemokrat._p']=round(PS_KOMPLET_17_96[,'PAR_2017_4_Českástr.sociálnědemokrat._p'],2)
PS_KOMPLET_17_96[,'PAR_2017_21_ANO2011_p']=round(PS_KOMPLET_17_96[,'PAR_2017_21_ANO2011_p'],2)
PS_KOMPLET_17_96[,'PAR10_9_Česká_strsociálně_demokrat_p']=round(PS_KOMPLET_17_96[,'PAR10_9_Česká_strsociálně_demokrat_p'],2)


#Standardizace vybraných proměnných - dataset FINAL_17
FINAL_17$sd_PAR_2017_29_OkamuraSPD_p <- scale(FINAL_17$PAR_2017_29_OkamuraSPD_p)
FINAL_17$sd_PAR_2017_21_ANO2011_p <- scale(FINAL_17$PAR_2017_21_ANO2011_p)
FINAL_17$sd_PAR_2017_8_KSCM_p <- scale(FINAL_17$PAR_2017_8_KSCM_p)
FINAL_17$sd_PAR_2017_25_ceskastrananarodnesocialni_p <- scale(FINAL_17$PAR_2017_25_ceskastrananarodnesocialni_p)

FINAL_17$sd_PAR13_21_KSCM_p <- scale(FINAL_17$PAR13_21_KSCM_p)
FINAL_17$sd_PAR13_20_ANO_2011_p <- scale(FINAL_17$PAR13_20_ANO_2011_p)
FINAL_17$sd_PAR13_17_Okamura_USVIT_p <- scale(FINAL_17$PAR13_17_Okamura_USVIT_p)
FINAL_17$sd_PAR13_1_CSSD_p <- scale(FINAL_17$PAR13_1_CSSD_p)

FINAL_17$sd_Vzd_zakl_A <- scale(FINAL_17$Vzd_zakl_A)
FINAL_17$sd_Eko_podnikatel_A <- scale(FINAL_17$Eko_podnikatel_A)
FINAL_17$sd_Eko_nezam_A <- scale(FINAL_17$Eko_nezam_A)
FINAL_17$sd_Rodst_rozv_A <- scale(FINAL_17$Rodst_rozv_A)
FINAL_17$sd_Nab_rimsko_A <- scale(FINAL_17$Nab_rimsko_A)
FINAL_17$sd_Vzd_vysok_A <- scale(FINAL_17$Vzd_vysok_A)
FINAL_17$sd_Eko_podnikatel_A <- scale(FINAL_17$Eko_podnikatel_A)
FINAL_17$sd_Eko_nezam_A <- scale(FINAL_17$Eko_nezam_A)
FINAL_17$sd_Nepracduch_A <- scale(FINAL_17$Nepracduch_A)
FINAL_17$sd_Vzd_zakl_A <- scale(FINAL_17$Vzd_zakl_A)
FINAL_17$sd_Sum_v0_14_A <- scale(FINAL_17$Sum_v0_14_A)
FINAL_17$sd_CSSD_ROZDIL_P <- scale(FINAL_17$CSSD_ROZDIL_P)
FINAL_17$sd_KDU_CSL_ROZDIL_P <- scale(FINAL_17$KDU_CSL_ROZDIL_P)
FINAL_17$sd_TOP09_ROZDIL_P <- scale(FINAL_17$TOP09_ROZDIL_P)
FINAL_17$sd_ANO_ROZDIL_P <- scale(FINAL_17$ANO_ROZDIL_P)
FINAL_17$sd_ODS_ROZDIL_P <- scale(FINAL_17$ODS_ROZDIL_P)
FINAL_17$sd_PIRATI_ROZDIL_P <- scale(FINAL_17$PIRATI_ROZDIL_P)
FINAL_17$sd_SPD_USVIT_ROZDIL_P <- scale(FINAL_17$SPD_USVIT_ROZDIL_P)
FINAL_17$sd_KSCM_ROZDIL_P <- scale(FINAL_17$KSCM_ROZDIL_P)

FINAL_17$sd_INHABITANTS_2018_LOGGED <- scale(FINAL_17$INHABITANTS_2018_LOGGED)
FINAL_17$sd_EXEKUCE_OBCE <- scale(FINAL_17$EXEKUCE_OBCE)



##############################################################################################################################
##############################################################################################################################
##############################################################################################################################


#Regrese ČSSD
MODEL_CSSD_17 <- lm(PAR_2017_4_Českástr.sociálnědemokrat._p ~ Vzd_vysok_A + Eko_podnikatel_A + Eko_nezam_A + Nepracduch_A + Nab_rimsko_A, data=PS_KOMPLET_17_96, weights=PAR_2017_VOL_SEZNAM)
summary(MODEL_CSSD_17)

MODEL_CSSD_13 <- lm(PAR13_1_Česká_strsociálně_demokrat_p ~ Vzd_vysok_A + Eko_podnikatel_A + Eko_nezam_A + Nepracduch_A + Nab_rimsko_A, data=PS_KOMPLET_17_96, weights=PAR_2017_VOL_SEZNAM)
summary(MODEL_CSSD_13) 

MODEL_CSSD_10 <- lm(PAR10_9_Česká_strsociálně_demokrat_p ~ Vzd_vysok_A + Eko_podnikatel_A + Eko_nezam_A + Nepracduch_A + Nab_rimsko_A, data=PS_KOMPLET_17_96, weights=PAR_2017_VOL_SEZNAM)
summary(MODEL_CSSD_10) 

MODEL_CSSD_06 <- lm(PAR06_10_Česká_strsociálně_demokrat_p ~ Vzd_vysok_A + Eko_podnikatel_A + Eko_nezam_A + Nepracduch_A + Nab_rimsko_A, data=PS_KOMPLET_17_96, weights=PAR_2017_VOL_SEZNAM)
summary(MODEL_CSSD_06) 

MODEL_CSSD_02 <- lm(PAR02_3_Česká_strsociálně_demokrat_p ~ Vzd_vysok_A + Eko_podnikatel_A + Eko_nezam_A + Nepracduch_A + Nab_rimsko_A, data=PS_KOMPLET_17_96, weights=PAR_2017_VOL_SEZNAM)
summary(MODEL_CSSD_02) 

MODEL_CSSD_98 <- lm(PAR98_7_Česká_strsociálně_demokrat_p ~ Vzd_vysok_A + Eko_podnikatel_A + Eko_nezam_A + Nepracduch_A + Nab_rimsko_A, data=PS_KOMPLET_17_96, weights=PAR_2017_VOL_SEZNAM)
summary(MODEL_CSSD_98) 

MODEL_CSSD_96 <- lm(PAR96_2_Česká_strsociáldemokratická_p ~ Vzd_vysok_A + Eko_podnikatel_A + Eko_nezam_A + Nepracduch_A + Nab_rimsko_A, data=PS_KOMPLET_17_96, weights=PAR_2017_VOL_SEZNAM)
summary(MODEL_CSSD_96) 

#Zápis do tabulky html 
library(texreg)
htmlreg(list(MODEL_CSSD_17, MODEL_CSSD_13, MODEL_CSSD_10, MODEL_CSSD_06, MODEL_CSSD_02, MODEL_CSSD_98, MODEL_CSSD_96), # názvy modelů
        file = "./MODELY_CSSD.html", # output file
        custom.model.names = c("ČSSD 2017" , "ČSSD 2013", "ČSSD 2010", "ČSSD 2006", "ČSSD 2002" , "ČSSD 1998", "ČSSD 1996"), # název modelů
        reorder.coef = c( 2, 3, 4, 5, 6,1),
        custom.coef.names = c("(Konstanta)","Vysokoškoláci" , "Podnikatelé", "Nezaměstnaní" , "Důchodci", "Katolící"), # názvy proměnných
        custom.gof.names = c("R2","Adj.R2","Počet pozorování", "RMSE"), # názvy proměnných
        digits = 3, # decimal places for coefficients and SEs
        single.row = FALSE) # display BETA and SEs under each other

#vytvoření grafu regresních koeficientů - nestandardizovaných
library(coefplot)
multiplot(MODEL_CSSD_17, MODEL_CSSD_13, MODEL_CSSD_10, MODEL_CSSD_06, MODEL_CSSD_02, MODEL_CSSD_98, MODEL_CSSD_96, intercept = FALSE,  title = "Volební podpora ČSSD", xlab = "Hodnota koeficientu - nestandardizovaná  + 95% k.i.",
          ylab = "Název  koeficientu", sort = c("magnitude"), pointSize = 3, color = "blue", cex = 5, 
          newNames=c(Vzd_vysok_A="Vysokoškoláci",
                     Eko_podnikatel_A="Podnikatelé",
                     Eko_nezam_A="Nezaměstnaní",
                     Nepracduch_A="Důchodci",
                     Nab_rimsko_A="Katolící"))

#Standardizované modely
SD_MODEL_CSSD_17 <- lm(sd_PAR_2017_4_Českástr.sociálnědemokrat._p ~ ZVzd_vysok_A + ZEko_podnikatel_A + ZEko_nezam_A + ZNepracduch_A + ZNab_rimsko_A + sd_EXEKUCE_OBCE + sd_INHABITANTS_2018_LOGGED, data=PS_KOMPLET_17_96, weights=PAR_2017_VOL_SEZNAM)
summary(SD_MODEL_CSSD_17)

SD_MODEL_CSSD_13 <- lm(sd_PAR13_1_Česká_strsociálně_demokrat_p ~ ZVzd_vysok_A + ZEko_podnikatel_A + ZEko_nezam_A + ZNepracduch_A + ZNab_rimsko_A + sd_EXEKUCE_OBCE + sd_INHABITANTS_2018_LOGGED, data=PS_KOMPLET_17_96, weights=PAR_2017_VOL_SEZNAM)
summary(SD_MODEL_CSSD_13) 

SD_MODEL_CSSD_10 <- lm(sd_PAR10_9_Česká_strsociálně_demokrat_p ~ ZVzd_vysok_A + ZEko_podnikatel_A + ZEko_nezam_A + ZNepracduch_A + ZNab_rimsko_A + sd_EXEKUCE_OBCE + sd_INHABITANTS_2018_LOGGED, data=PS_KOMPLET_17_96, weights=PAR_2017_VOL_SEZNAM)
summary(SD_MODEL_CSSD_10) 

SD_MODEL_CSSD_06 <- lm(sd_PAR06_10_Česká_strsociálně_demokrat_p ~ ZVzd_vysok_A + ZEko_podnikatel_A + ZEko_nezam_A + ZNepracduch_A + ZNab_rimsko_A + sd_EXEKUCE_OBCE + sd_INHABITANTS_2018_LOGGED, data=PS_KOMPLET_17_96, weights=PAR_2017_VOL_SEZNAM)
summary(SD_MODEL_CSSD_06) 

SD_MODEL_CSSD_02 <- lm(sd_PAR02_3_Česká_strsociálně_demokrat_p ~ ZVzd_vysok_A + ZEko_podnikatel_A + ZEko_nezam_A + ZNepracduch_A + ZNab_rimsko_A + sd_EXEKUCE_OBCE + sd_INHABITANTS_2018_LOGGED, data=PS_KOMPLET_17_96, weights=PAR_2017_VOL_SEZNAM)
summary(SD_MODEL_CSSD_02) 

SD_MODEL_CSSD_98 <- lm(sd_PAR98_7_Česká_strsociálně_demokrat_p ~ ZVzd_vysok_A + ZEko_podnikatel_A + ZEko_nezam_A + ZNepracduch_A + ZNab_rimsko_A + sd_EXEKUCE_OBCE + sd_INHABITANTS_2018_LOGGED, data=PS_KOMPLET_17_96, weights=PAR_2017_VOL_SEZNAM)
summary(SD_MODEL_CSSD_98) 

SD_MODEL_CSSD_96 <- lm(sd_PAR96_2_Česká_strsociáldemokratická_p ~ ZVzd_vysok_A + ZEko_podnikatel_A + ZEko_nezam_A + ZNepracduch_A + ZNab_rimsko_A + sd_EXEKUCE_OBCE + sd_INHABITANTS_2018_LOGGED, data=PS_KOMPLET_17_96, weights=PAR_2017_VOL_SEZNAM)
summary(SD_MODEL_CSSD_96) 


#Zápis do tabulky html 
library(texreg)
htmlreg(list(SD_MODEL_CSSD_17, SD_MODEL_CSSD_13, SD_MODEL_CSSD_10, SD_MODEL_CSSD_06, SD_MODEL_CSSD_02,
             SD_MODEL_CSSD_98, SD_MODEL_CSSD_96), # názvy modelů
        file = "./MODELY_CSSD_SD.html", # output file
        custom.model.names = c("ČSSD 2017" , "ČSSD 2013", "ČSSD 2010", "ČSSD 2006", "ČSSD 2002" , "ČSSD 1998", "ČSSD 1996"), # název modelů
        reorder.coef = c( 2, 3, 4, 5, 6,1),
        custom.coef.names = c("(Konstanta)","Vysokoškoláci" , "Podnikatelé", "Nezaměstnaní" , "Důchodci", "Katolící"), # názvy proměnných
        custom.gof.names = c("R2","Adj.R2","Počet pozorování", "RMSE"), # názvy proměnných
        digits = 3, # decimal places for coefficients and SEs
        single.row = FALSE) # display BETA and SEs under each other

#vytvoření grafu regresních koeficientů - standardizovaných
library(coefplot)
multiplot(SD_MODEL_CSSD_17, SD_MODEL_CSSD_13, SD_MODEL_CSSD_10, SD_MODEL_CSSD_06, SD_MODEL_CSSD_02,
          SD_MODEL_CSSD_98, SD_MODEL_CSSD_96, intercept = FALSE,  title = "Volební podpora ČSSD", xlab = "Hodnota koeficientu - standardizovaná  + 95% k.i.",
          ylab = "Název  koeficientu", sort = c("magnitude"), pointSize = 3, color = "blue", cex = 5, 
          newNames=c(ZVzd_vysok_A="Vysokoškoláci",
                     ZEko_podnikatel_A="Podnikatelé",
                     ZEko_nezam_A="Nezaměstnaní",
                     ZNepracduch_A="Důchodci",
                     ZNab_rimsko_A="Katolící"))

#vytvoření grafu regresních koeficientů - standardizovaných - jen od roku 2006
library(coefplot)
multiplot(SD_MODEL_CSSD_17, SD_MODEL_CSSD_13, SD_MODEL_CSSD_10, SD_MODEL_CSSD_06, intercept = FALSE,  title = "Volební podpora ČSSD", xlab = "Hodnota koeficientu - standardizovaná  + 95% k.i.",
          ylab = "Název  koeficientu", sort = c("magnitude"), pointSize = 3, color = "blue", cex = 5, 
          newNames=c(ZVzd_vysok_A="Vysokoškoláci",
                     ZEko_podnikatel_A="Podnikatelé",
                     ZEko_nezam_A="Nezaměstnaní",
                     ZNepracduch_A="Důchodci"))

#vytvoření grafu regresních koeficientů - standardizovaných - jen od roku 2013
SD_MODEL_CSSD_17 <- lm(sd_PAR_2017_4_Českástr.sociálnědemokrat._p ~ ZVzd_vysok_A + ZEko_podnikatel_A + ZEko_nezam_A + ZNepracduch_A, data=PS_KOMPLET_17_96, weights=PAR_2017_VOL_SEZNAM)
summary(SD_MODEL_CSSD_17)

SD_MODEL_CSSD_13 <- lm(sd_PAR13_1_Česká_strsociálně_demokrat_p ~ ZVzd_vysok_A + ZEko_podnikatel_A + ZEko_nezam_A + ZNepracduch_A, data=PS_KOMPLET_17_96, weights=PAR_2017_VOL_SEZNAM)
summary(SD_MODEL_CSSD_13) 


# v článku!
library(coefplot)
multiplot(SD_MODEL_CSSD_17, SD_MODEL_CSSD_13, intercept = FALSE,  title = "DV: ČSSD - vote share in %", xlab = "Standardized coefficients  + 95% c.i.",
          ylab = "Name", sort = c("magnitude"), pointSize = 3, color = "blue", cex = 5,  names=c("2017",  "2013"),
          newNames=c(ZVzd_vysok_A="College degree (%)",
                     ZEko_podnikatel_A="Self-employed (%)",
                     ZEko_nezam_A="Unemployed (%)",
                     ZNepracduch_A="65+ years (%)")) +
  geom_vline(xintercept=0, linetype="dashed", color = "red", size=1.1) + theme(legend.position="bottom") + theme(legend.position="bottom")  + theme(plot.title = element_text(size=12,face="bold")) +
  theme(axis.text=element_text(size=13), axis.title=element_text(size=10,face="bold")) + scale_color_manual(values=c("black","blue"))





#Analýza ztrát ČSSD
ROZ_CSSD_98_96 <- lm(ROZ_CSSD_98_96 ~ Vzd_vysok_A + Eko_podnikatel_A + Eko_nezam_A + Nepracduch_A + Nab_rimsko_A, data=PS_KOMPLET_17_96, weights=PAR_2017_VOL_SEZNAM)
summary(ROZ_CSSD_98_96) 

ROZ_CSSD_02_98 <- lm(ROZ_CSSD_02_98 ~ Vzd_vysok_A + Eko_podnikatel_A + Eko_nezam_A + Nepracduch_A + Nab_rimsko_A, data=PS_KOMPLET_17_96, weights=PAR_2017_VOL_SEZNAM)
summary(ROZ_CSSD_02_98) 

ROZ_CSSD_06_02 <- lm(ROZ_CSSD_06_02 ~ Vzd_vysok_A + Eko_podnikatel_A + Eko_nezam_A + Nepracduch_A + Nab_rimsko_A, data=PS_KOMPLET_17_96, weights=PAR_2017_VOL_SEZNAM)
summary(ROZ_CSSD_06_02) 

ROZ_CSSD_10_06 <- lm(ROZ_CSSD_10_06 ~ Vzd_vysok_A + Eko_podnikatel_A + Eko_nezam_A + Nepracduch_A + Nab_rimsko_A, data=PS_KOMPLET_17_96, weights=PAR_2017_VOL_SEZNAM)
summary(ROZ_CSSD_10_06) 

ROZ_CSSD_13_10 <- lm(ROZ_CSSD_13_10 ~ Vzd_vysok_A + Eko_podnikatel_A + Eko_nezam_A + Nepracduch_A + Nab_rimsko_A, data=PS_KOMPLET_17_96, weights=PAR_2017_VOL_SEZNAM)
summary(ROZ_CSSD_13_10) 

ROZ_CSSD_17_13 <- lm(ROZ_CSSD_17_13 ~ Vzd_vysok_A + Eko_podnikatel_A + Eko_nezam_A + Nepracduch_A + Nab_rimsko_A, data=PS_KOMPLET_17_96, weights=PAR_2017_VOL_SEZNAM)
summary(ROZ_CSSD_17_13) 


#Zápis do tabulky html 
library(texreg)
htmlreg(list(ROZ_CSSD_98_96, ROZ_CSSD_02_98, ROZ_CSSD_06_02, ROZ_CSSD_10_06, ROZ_CSSD_13_10,
             ROZ_CSSD_17_13), # názvy modelů
        file = "./ROZ_MODELY_CSSD.html", # output file
        custom.model.names = c("ČSSD 1998-96", "ČSSD 2002-98",  "ČSSD 2006-02",  "ČSSD 2010-06",  "ČSSD 2013-10", "ČSSD 2017-13"), # název modelů
        reorder.coef = c( 2, 3, 4, 5, 6,1),
        custom.coef.names = c("(Konstanta)","Vysokoškoláci" , "Podnikatelé", "Nezaměstnaní" , "Důchodci", "Katolící"), # názvy proměnných
        custom.gof.names = c("R2","Adj.R2","Počet pozorování", "RMSE"), # názvy proměnných
        digits = 3, # decimal places for coefficients and SEs
        single.row = FALSE) # display BETA and SEs under each other


#Rozdíl u vybraných voleb s volbami 2017
ROZ_CSSD_17_06 <- lm(ROZ_CSSD_17_06 ~ Vzd_vysok_A + Eko_podnikatel_A + Eko_nezam_A + Nepracduch_A + Nab_rimsko_A, data=PS_KOMPLET_17_96, weights=PAR_2017_VOL_SEZNAM)
summary(ROZ_CSSD_17_06) 
ROZ_CSSD_17_96 <- lm(ROZ_CSSD_17_06 ~ Vzd_vysok_A + Eko_podnikatel_A + Eko_nezam_A + Nepracduch_A + Nab_rimsko_A, data=PS_KOMPLET_17_96, weights=PAR_2017_VOL_SEZNAM)
summary(ROZ_CSSD_17_06) 

#Analýza KSČM

#Regrese KSČM
MODEL_KSCM_17 <- lm(PAR_2017_8_Komunistickástr.ČechaMoravy_p ~ Vzd_vysok_A + Eko_podnikatel_A + Eko_nezam_A + Nepracduch_A + Nab_rimsko_A, data=PS_KOMPLET_17_96, weights=PAR_2017_VOL_SEZNAM)
summary(MODEL_KSCM_17)

MODEL_KSCM_13 <- lm(PAR13_21_Komunistická_strČech_a_Moravy_p ~ Vzd_vysok_A + Eko_podnikatel_A + Eko_nezam_A + Nepracduch_A + Nab_rimsko_A, data=PS_KOMPLET_17_96, weights=PAR_2017_VOL_SEZNAM)
summary(MODEL_KSCM_13) 

MODEL_KSCM_10 <- lm(PAR10_6_Komunistická_strČech_a_Moravy_p ~ Vzd_vysok_A + Eko_podnikatel_A + Eko_nezam_A + Nepracduch_A + Nab_rimsko_A, data=PS_KOMPLET_17_96, weights=PAR_2017_VOL_SEZNAM)
summary(MODEL_KSCM_10) 

MODEL_KSCM_06 <- lm(PAR06_20_Komunistická_strČech_a_Moravy_p ~ Vzd_vysok_A + Eko_podnikatel_A + Eko_nezam_A + Nepracduch_A + Nab_rimsko_A, data=PS_KOMPLET_17_96, weights=PAR_2017_VOL_SEZNAM)
summary(MODEL_KSCM_06) 

MODEL_KSCM_02 <- lm(PAR02_23_Komunistická_strČech_a_Moravy_p ~ Vzd_vysok_A + Eko_podnikatel_A + Eko_nezam_A + Nepracduch_A + Nab_rimsko_A, data=PS_KOMPLET_17_96, weights=PAR_2017_VOL_SEZNAM)
summary(MODEL_KSCM_02) 

MODEL_KSCM_98 <- lm(PAR98_9_Komunistická_strČech_a_Moravy_p ~ Vzd_vysok_A + Eko_podnikatel_A + Eko_nezam_A + Nepracduch_A + Nab_rimsko_A, data=PS_KOMPLET_17_96, weights=PAR_2017_VOL_SEZNAM)
summary(MODEL_KSCM_98) 

MODEL_KSCM_96 <- lm(PAR96_15_Komunistická_strČech_a_Moravy_p ~ Vzd_vysok_A + Eko_podnikatel_A + Eko_nezam_A + Nepracduch_A + Nab_rimsko_A, data=PS_KOMPLET_17_96, weights=PAR_2017_VOL_SEZNAM)
summary(MODEL_KSCM_96) 

library(texreg)
htmlreg(list(MODEL_KSCM_17, MODEL_KSCM_13, MODEL_KSCM_10, MODEL_KSCM_06, MODEL_KSCM_02, MODEL_KSCM_98, MODEL_KSCM_96), # názvy modelů
        file = "./MODELY_KSCM.html", # output file
        custom.model.names = c("KSČM 2017" , "KSČM 2013", "KSČM 2010", "KSČM 2006", "KSČM 2002" , "KSČM 1998", "KSČM 1996"), # název modelů
        reorder.coef = c( 2, 3, 4, 5, 6,1),
        custom.coef.names = c("(Konstanta)","Vysokoškoláci" , "Podnikatelé", "Nezaměstnaní" , "Důchodci", "Katolící"), # názvy proměnných
        custom.gof.names = c("R2","Adj.R2","Počet pozorování", "RMSE"), # názvy proměnných
        digits = 3, # decimal places for coefficients and SEs
        single.row = FALSE) # display BETA and SEs under each other

#vytvoření grafu regresních koeficientů - nestandardizovaných
library(coefplot)
multiplot(MODEL_KSCM_17, MODEL_KSCM_13, KSCM_CSSD_10, MODEL_KSCM_06, MODEL_KSCM_02, MODEL_KSCM_98, MODEL_KSCM_96, intercept = FALSE,  title = "Volební podpora ČSSD", xlab = "Hodnota koeficientu - nestandardizovaná  + 95% k.i.",
          ylab = "Název  koeficientu", sort = c("magnitude"), pointSize = 3, color = "blue", cex = 5, 
          newNames=c(Vzd_vysok_A="Vysokoškoláci",
                     Eko_podnikatel_A="Podnikatelé",
                     Eko_nezam_A="Nezaměstnaní",
                     Nepracduch_A="Důchodci",
                     Nab_rimsko_A="Katolící"))

# Standardizované modely

SD_MODEL_KSCM_17 <- lm(sd_PAR_2017_8_Komunistickástr.ČechaMoravy_p ~ Vzd_vysok_A + Eko_podnikatel_A + Eko_nezam_A + Nepracduch_A + Nab_rimsko_A, data=PS_KOMPLET_17_96, weights=PAR_2017_VOL_SEZNAM)
summary(SD_MODEL_KSCM_17)

SD_MODEL_KSCM_13 <- lm(sd_PAR13_21_Komunistická_strČech_a_Moravy_p ~ Vzd_vysok_A + Eko_podnikatel_A + Eko_nezam_A + Nepracduch_A + Nab_rimsko_A, data=PS_KOMPLET_17_96, weights=PAR_2017_VOL_SEZNAM)
summary(SD_MODEL_KSCM_13) 

SD_MODEL_KSCM_10 <- lm(sd_PAR10_6_Komunistická_strČech_a_Moravy_p ~ Vzd_vysok_A + Eko_podnikatel_A + Eko_nezam_A + Nepracduch_A + Nab_rimsko_A, data=PS_KOMPLET_17_96, weights=PAR_2017_VOL_SEZNAM)
summary(SD_MODEL_KSCM_10) 

SD_MODEL_KSCM_06 <- lm(sd_PAR06_20_Komunistická_strČech_a_Moravy_p ~ Vzd_vysok_A + Eko_podnikatel_A + Eko_nezam_A + Nepracduch_A + Nab_rimsko_A, data=PS_KOMPLET_17_96, weights=PAR_2017_VOL_SEZNAM)
summary(SD_MODEL_KSCM_06) 

SD_MODEL_KSCM_02 <- lm(sd_PAR02_23_Komunistická_strČech_a_Moravy_p ~ Vzd_vysok_A + Eko_podnikatel_A + Eko_nezam_A + Nepracduch_A + Nab_rimsko_A, data=PS_KOMPLET_17_96, weights=PAR_2017_VOL_SEZNAM)
summary(SD_MODEL_KSCM_02) 

SD_MODEL_KSCM_98 <- lm(sd_PAR98_9_Komunistická_strČech_a_Moravy_p ~ Vzd_vysok_A + Eko_podnikatel_A + Eko_nezam_A + Nepracduch_A + Nab_rimsko_A, data=PS_KOMPLET_17_96, weights=PAR_2017_VOL_SEZNAM)
summary(SD_MODEL_KSCM_98) 

SD_MODEL_KSCM_96 <- lm(sd_PAR96_15_Komunistická_strČech_a_Moravy_p ~ Vzd_vysok_A + Eko_podnikatel_A + Eko_nezam_A + Nepracduch_A + Nab_rimsko_A, data=PS_KOMPLET_17_96, weights=PAR_2017_VOL_SEZNAM)
summary(SD_MODEL_KSCM_96) 

library(texreg)
htmlreg(list(SD_MODEL_KSCM_17, SD_MODEL_KSCM_13, SD_MODEL_KSCM_10, SD_MODEL_KSCM_06, SD_MODEL_KSCM_02, SD_MODEL_KSCM_98, SD_MODEL_KSCM_96), # názvy modelů
        file = "./MODELY_KSCM_SD.html", # output file
        custom.model.names = c("KSČM 2017" , "KSČM 2013", "KSČM 2010", "KSČM 2006", "KSČM 2002" , "KSČM 1998", "KSČM 1996"), # název modelů
        reorder.coef = c( 2, 3, 4, 5, 6,1),
        custom.coef.names = c("(Konstanta)","Vysokoškoláci" , "Podnikatelé", "Nezaměstnaní" , "Důchodci", "Katolící"), # názvy proměnných
        custom.gof.names = c("R2","Adj.R2","Počet pozorování", "RMSE"), # názvy proměnných
        digits = 3, # decimal places for coefficients and SEs
        single.row = FALSE) # display BETA and SEs under each other

#vytvoření grafu regresních koeficientů - nestandardizovaných
library(coefplot)
multiplot(SD_MODEL_KSCM_17, SD_MODEL_KSCM_13, SD_KSCM_CSSD_10, SD_MODEL_KSCM_06, SD_MODEL_KSCM_02, SD_MODEL_KSCM_98, SD_MODEL_KSCM_96, intercept = FALSE,  title = "Volební podpora ČSSD", xlab = "Hodnota koeficientu - nestandardizovaná  + 95% k.i.",
          ylab = "Název  koeficientu", sort = c("magnitude"), pointSize = 3, color = "blue", cex = 5, 
          newNames=c(Vzd_vysok_A="Vysokoškoláci",
                     Eko_podnikatel_A="Podnikatelé",
                     Eko_nezam_A="Nezaměstnaní",
                     Nepracduch_A="Důchodci",
                     Nab_rimsko_A="Katolící"))




#Analýza ztrát KSČM
ROZ_KSCM_98_96 <- lm(ROZ_KSCM_98_96 ~ Vzd_vysok_A + Eko_podnikatel_A + Eko_nezam_A + Nepracduch_A + Nab_rimsko_A, data=PS_KOMPLET_17_96, weights=PAR_2017_VOL_SEZNAM)
summary(ROZ_KSCM_98_96) 

ROZ_KSCM_02_98 <- lm(ROZ_KSCM_02_98 ~ Vzd_vysok_A + Eko_podnikatel_A + Eko_nezam_A + Nepracduch_A + Nab_rimsko_A, data=PS_KOMPLET_17_96, weights=PAR_2017_VOL_SEZNAM)
summary(ROZ_KSCM_02_98) 

ROZ_KSCM_06_02 <- lm(ROZ_KSCM_06_02 ~ Vzd_vysok_A + Eko_podnikatel_A + Eko_nezam_A + Nepracduch_A + Nab_rimsko_A, data=PS_KOMPLET_17_96, weights=PAR_2017_VOL_SEZNAM)
summary(ROZ_KSCM_06_02) 

ROZ_KSCM_10_06 <- lm(ROZ_KSCM_10_06 ~ Vzd_vysok_A + Eko_podnikatel_A + Eko_nezam_A + Nepracduch_A + Nab_rimsko_A, data=PS_KOMPLET_17_96, weights=PAR_2017_VOL_SEZNAM)
summary(ROZ_KSCM_10_06) 

ROZ_KSCM_13_10 <- lm(ROZ_KSCM_13_10 ~ Vzd_vysok_A + Eko_podnikatel_A + Eko_nezam_A + Nepracduch_A + Nab_rimsko_A, data=PS_KOMPLET_17_96, weights=PAR_2017_VOL_SEZNAM)
summary(ROZ_KSCM_13_10) 

ROZ_KSCM_17_13 <- lm(ROZ_KSCM_17_13 ~ Vzd_vysok_A + Eko_podnikatel_A + Eko_nezam_A + Nepracduch_A + Nab_rimsko_A, data=PS_KOMPLET_17_96, weights=PAR_2017_VOL_SEZNAM)
summary(ROZ_KSCM_17_13) 

#Zápis do tabulky html 
library(texreg)
htmlreg(list(ROZ_KSCM_98_96, ROZ_KSCM_02_98, ROZ_KSCM_06_02, ROZ_KSCM_10_06, ROZ_KSCM_13_10,
             ROZ_KSCM_17_13), # názvy modelů
        file = "./ROZ_MODELY_KSČM.html", # output file
        custom.model.names = c("KSČM 1998-96", "KSČM 2002-98",  "KSČM 2006-02",  "KSČM 2010-06",  "KSČM 2013-10", "KSČM 2017-13"), # název modelů
        reorder.coef = c( 2, 3, 4, 5, 6,1),
        custom.coef.names = c("(Konstanta)","Vysokoškoláci" , "Podnikatelé", "Nezaměstnaní" , "Důchodci", "Katolící"), # názvy proměnných
        custom.gof.names = c("R2","Adj.R2","Počet pozorování", "RMSE"), # názvy proměnných
        digits = 3, # decimal places for coefficients and SEs
        single.row = FALSE) # display BETA and SEs under each other



# Analýza ANO a SPD + doplnění komplet dat za KSČM a ČSSD - nestandardizované


ANO_17 <- lm(PAR_2017_21_ANO2011_p ~ Vzd_vysok_A + Eko_podnikatel_A + Eko_nezam_A + Nepracduch_A + Nab_rimsko_A + EXEKUCE_OBCE + INHABITANTS_2018_LOGGED, data=FINAL_17, weights=VOL_SEZNAM_17)
summary(ANO_17)

SPD_USVIT_17 <- lm(PAR_2017_29_OkamuraSPD_p ~ Vzd_vysok_A + Eko_podnikatel_A + Eko_nezam_A + Nepracduch_A + Nab_rimsko_A + EXEKUCE_OBCE + INHABITANTS_2018_LOGGED, data=FINAL_17, weights=VOL_SEZNAM_17)
summary(SPD_USVIT_17)

CSSD_17 <- lm(PAR_2017_25_ceskastrananarodnesocialni_p ~ Vzd_vysok_A + Eko_podnikatel_A + Eko_nezam_A + Nepracduch_A + Nab_rimsko_A + EXEKUCE_OBCE + INHABITANTS_2018_LOGGED, data=FINAL_17, weights=VOL_SEZNAM_17)
summary(CSSD_17)

KSCM_17 <- lm(PAR_2017_8_KSCM_p ~ Vzd_vysok_A + Eko_podnikatel_A + Eko_nezam_A + Nepracduch_A + Nab_rimsko_A + EXEKUCE_OBCE + INHABITANTS_2018_LOGGED, data=FINAL_17, weights=VOL_SEZNAM_17)
summary(KSCM_17)

ANO_13 <- lm(PAR13_20_ANO_2011_p ~ Vzd_vysok_A + Eko_podnikatel_A + Eko_nezam_A + Nepracduch_A + Nab_rimsko_A + EXEKUCE_OBCE + INHABITANTS_2018_LOGGED, data=FINAL_17, weights=VOL_SEZNAM_17)
summary(ANO_13)

SPD_USVIT_13 <- lm(PAR13_17_Okamura_USVIT_p ~ Vzd_vysok_A + Eko_podnikatel_A + Eko_nezam_A + Nepracduch_A + Nab_rimsko_A + EXEKUCE_OBCE + INHABITANTS_2018_LOGGED, data=FINAL_17, weights=VOL_SEZNAM_17)
summary(SPD_USVIT_13)

CSSD_13 <- lm(PAR13_1_CSSD_p ~ Vzd_vysok_A + Eko_podnikatel_A + Eko_nezam_A + Nepracduch_A + Nab_rimsko_A + EXEKUCE_OBCE + INHABITANTS_2018_LOGGED, data=FINAL_17, weights=VOL_SEZNAM_17)
summary(CSSD_17)

KSCM_13 <- lm(PAR13_21_KSCM_p ~ Vzd_vysok_A + Eko_podnikatel_A + Eko_nezam_A + Nepracduch_A + Nab_rimsko_A + EXEKUCE_OBCE + INHABITANTS_2018_LOGGED, data=FINAL_17, weights=VOL_SEZNAM_17)
summary(KSCM_13)



#Zápis do tabulky html  -  TABULKA V ČLÁNKU!!!!!!!!!!!!
library(texreg)
htmlreg(list(ANO_17, SPD_USVIT_17, CSSD_17, KSCM_17, ANO_13, SPD_USVIT_13, CSSD_13, KSCM_13), # názvy modelů
        file = "./STRANY_2017_2013.html", # output file
        custom.model.names = c("ANO 2017", "SPD 2017", "CSSD 2017", "KSČM 2017",
                               "ANO 2013", "SPD 2013", "CSSD 2013", "KSČM 2013"), # název modelů
        reorder.coef = c( 2, 3, 4, 5, 6, 7, 8, 1),
        custom.coef.names = c("(Constant)","College degree (%)" , "Self-employed (%)", "Unemployed (%)" , "65+ years (%)", "Catholics (%)",
                              "Distraints (%)", "Number of inhabitants (log)"), # názvy proměnných
        custom.gof.names = c("R2","Adj.R2","Number of cases", "RMSE"), # názvy proměnných
        digits = 3, # decimal places for coefficients and SEs
        single.row = FALSE) # display BETA and SEs under each other


sd_ANO_17 <- lm(sd_PAR_2017_21_ANO2011_p ~ sd_Vzd_vysok_A + sd_Eko_podnikatel_A + sd_Eko_nezam_A + sd_Nepracduch_A + sd_Nab_rimsko_A + sd_EXEKUCE_OBCE + sd_INHABITANTS_2018_LOGGED, data=FINAL_17, weights=VOL_SEZNAM_17)
summary(sd_ANO_17)

sd_SPD_USVIT_17 <- lm(sd_PAR_2017_29_OkamuraSPD_p ~ sd_Vzd_vysok_A + sd_Eko_podnikatel_A + sd_Eko_nezam_A + sd_Nepracduch_A + sd_Nab_rimsko_A  + sd_EXEKUCE_OBCE + sd_INHABITANTS_2018_LOGGED, data=FINAL_17, weights=VOL_SEZNAM_17)
summary(sd_SPD_USVIT_17)

sd_CSSD_17 <- lm(sd_PAR_2017_25_ceskastrananarodnesocialni_p ~ sd_Vzd_vysok_A + sd_Eko_podnikatel_A + sd_Eko_nezam_A + sd_Nepracduch_A + sd_Nab_rimsko_A + sd_EXEKUCE_OBCE + sd_INHABITANTS_2018_LOGGED, data=FINAL_17, weights=VOL_SEZNAM_17)
summary(sd_CSSD_17)

sd_KSCM_17 <- lm(sd_PAR_2017_8_KSCM_p ~ sd_Vzd_vysok_A + sd_Eko_podnikatel_A + sd_Eko_nezam_A + sd_Nepracduch_A + sd_Nab_rimsko_A + sd_EXEKUCE_OBCE + sd_INHABITANTS_2018_LOGGED, data=FINAL_17, weights=VOL_SEZNAM_17)
summary(sd_KSCM_17)


sd_ANO_13 <- lm(sd_PAR13_20_ANO_2011_p ~  sd_Vzd_vysok_A + sd_Eko_podnikatel_A + sd_Eko_nezam_A + sd_Nepracduch_A + sd_Nab_rimsko_A + sd_EXEKUCE_OBCE + sd_INHABITANTS_2018_LOGGED, data=FINAL_17, weights=VOL_SEZNAM_17)
summary(sd_ANO_13)

sd_SPD_USVIT_13 <- lm(sd_PAR13_17_Okamura_USVIT_p ~  sd_Vzd_vysok_A + sd_Eko_podnikatel_A + sd_Eko_nezam_A + sd_Nepracduch_A + sd_Nab_rimsko_A + sd_EXEKUCE_OBCE + sd_INHABITANTS_2018_LOGGED, data=FINAL_17, weights=VOL_SEZNAM_17)
summary(sd_SPD_USVIT_13)

sd_CSSD_13 <- lm(sd_PAR13_1_CSSD_p ~  sd_Vzd_vysok_A + sd_Eko_podnikatel_A + sd_Eko_nezam_A + sd_Nepracduch_A + sd_Nab_rimsko_A + sd_EXEKUCE_OBCE + sd_INHABITANTS_2018_LOGGED, data=FINAL_17, weights=VOL_SEZNAM_17)
summary(sd_CSSD_17)

sd_KSCM_13 <- lm(sd_PAR13_21_KSCM_p ~  sd_Vzd_vysok_A + sd_Eko_podnikatel_A + sd_Eko_nezam_A + sd_Nepracduch_A + sd_Nab_rimsko_A + sd_EXEKUCE_OBCE + sd_INHABITANTS_2018_LOGGED, data=FINAL_17, weights=VOL_SEZNAM_17)
summary(sd_KSCM_13)

library(texreg)
htmlreg(list(sd_ANO_17, sd_SPD_USVIT_17, sd_CSSD_17, sd_KSCM_17, sd_ANO_13, sd_SPD_USVIT_13, sd_CSSD_13, sd_KSCM_13), # názvy modelů
        file = "./STRANY_2017_2013_SD.html", # output file
        custom.model.names = c("ANO 2017", "SPD 2017", "CSSD 2017", "KSČM 2017",
                               "ANO 2013", "SPD 2013", "CSSD 2013", "KSČM 2013"), # název modelů
        reorder.coef = c( 2, 3, 4, 5, 6, 7, 8, 1),
        custom.coef.names = c("(Konstanta)","Vysokoškoláci" , "Podnikatelé", "Nezaměstnaní" , "Důchodci", "Katolící",
                              "Distraints (%)", "Number of inhabitants (log)"), # názvy proměnných
        custom.gof.names = c("R2","Adj.R2","Počet pozorování", "RMSE"), # názvy proměnných
        digits = 3, # decimal places for coefficients and SEs
        single.row = FALSE) # display BETA and SEs under each other











# ROZDÍL

ROZ_ANO_17_13 <- lm(ANO_ROZDIL_P ~ Vzd_vysok_A + Eko_podnikatel_A + Eko_nezam_A + Nepracduch_A + Nab_rimsko_A + EXEKUCE_OBCE + INHABITANTS_2018_LOGGED, data=FINAL_17, weights=VOL_SEZNAM_17)
summary(ROZ_ANO_17_13)

ROZ_SPD_USVIT_17_13 <- lm(SPD_USVIT_ROZDIL_P ~ Vzd_vysok_A + Eko_podnikatel_A + Eko_nezam_A + Nepracduch_A + Nab_rimsko_A + EXEKUCE_OBCE + INHABITANTS_2018_LOGGED, data=FINAL_17, weights=VOL_SEZNAM_17)
summary(ROZ_SPD_USVIT_17_13)

ROZ_CSSD_17_13 <- lm(CSSD_ROZDIL_P ~ Vzd_vysok_A + Eko_podnikatel_A + Eko_nezam_A + Nepracduch_A + Nab_rimsko_A + EXEKUCE_OBCE + INHABITANTS_2018_LOGGED, data=FINAL_17, weights=VOL_SEZNAM_17)
summary(ROZ_CSSD_17_13)

ROZ_KSCM_17_13 <- lm(KSCM_ROZDIL_P ~ Vzd_vysok_A + Eko_podnikatel_A + Eko_nezam_A + Nepracduch_A + Nab_rimsko_A + EXEKUCE_OBCE + INHABITANTS_2018_LOGGED, data=FINAL_17, weights=VOL_SEZNAM_17)
summary(ROZ_KSCM_17_13)

#Zápis do tabulky html  -  TABULKA V ČLÁNKU!!!!!!!!!!!!
library(texreg)
htmlreg(list(ROZ_ANO_17_13, ROZ_SPD_USVIT_17_13, ROZ_CSSD_17_13, ROZ_KSCM_17_13), # názvy modelů
        file = "./ROZ_STRANY_2017_2013.html", # output file
        custom.model.names = c("ANO 2017-13", "SPD 2017-13", "CSSD 2017-13", "KSČM 2017-13"), # název modelů
        reorder.coef = c( 2, 3, 4, 5, 6, 7, 8, 1),
        custom.coef.names = c("(Constant)","College degree (%)" , "Self-employed (%)", "Unemployed (%)" , "65+ years (%)", "Catholics (%)",
                              "Distraints (%)", "Number of inhabitants (log)"), # názvy proměnných
        custom.gof.names = c("R2","Adj.R2","Number of cases", "RMSE"), # názvy proměnných
        digits = 3, # decimal places for coefficients and SEs
        single.row = FALSE) # display BETA and SEs under each other


# Analýza ANO a SPD + doplnění komplet dat za KSČM a ČSSD - nestandardizované a HLM - okres
ROZ_ANO_17_13 <- lmer(ANO_ROZDIL_P ~ Vzd_vysok_A + Eko_podnikatel_A + Eko_nezam_A + Nepracduch_A + Nab_rimsko_A + EXEKUCE_OBCE + INHABITANTS_2018_LOGGED + (1| OKRES), data=FINAL_17, weights=VOL_SEZNAM_17)
summary(ROZ_ANO_17_13)

ROZ_SPD_USVIT_17_13 <- lmer(SPD_USVIT_ROZDIL_P ~ Vzd_vysok_A + Eko_podnikatel_A + Eko_nezam_A + Nepracduch_A + Nab_rimsko_A + EXEKUCE_OBCE + INHABITANTS_2018_LOGGED + (1 + INHABITANTS_2018_LOGGED | KRAJ), data=FINAL_17, weights=VOL_SEZNAM_17)
summary(ROZ_SPD_USVIT_17_13)
ranef(ROZ_SPD_USVIT_17_13)

ROZ_CSSD_17_13 <- lmer(CSSD_ROZDIL_P ~ Vzd_vysok_A + Eko_podnikatel_A + Eko_nezam_A + Nepracduch_A + Nab_rimsko_A + EXEKUCE_OBCE + INHABITANTS_2018_LOGGED + (1| OKRES), data=FINAL_17, weights=VOL_SEZNAM_17)
summary(ROZ_CSSD_17_13)

ROZ_KSCM_17_13 <- lmer(KSCM_ROZDIL_P ~ Vzd_vysok_A + Eko_podnikatel_A + Eko_nezam_A + Nepracduch_A + Nab_rimsko_A + EXEKUCE_OBCE + INHABITANTS_2018_LOGGED + (1| OKRES), data=FINAL_17, weights=VOL_SEZNAM_17)
summary(ROZ_KSCM_17_13)

#nulové modely
NULL_ROZ_ANO_17_13 <- lmer(ANO_ROZDIL_P ~ 1 + (1| OKRES), data=FINAL_17)
summary(NULL_ROZ_ANO_17_13)
NULL_ROZ_SPD_USVIT_17_13 <- lmer(SPD_USVIT_ROZDIL_P ~ 1 + (1 | OKRES), data=FINAL_17)
summary(NULL_ROZ_SPD_USVIT_17_13)
NULL_ROZ_CSSD_17_13 <- lmer(CSSD_ROZDIL_P ~ 1 + (1 | OKRES), data=FINAL_17)
summary(NULL_ROZ_CSSD_17_13)
NULL_KSCM_ROZDIL_17_13 <- lmer(KSCM_ROZDIL_P ~ 1 + (1 | OKRES), data=FINAL_17)
summary(NULL_KSCM_ROZDIL_17_13)

#Zápis do tabulky html  -  
library(texreg)
htmlreg(list(ROZ_ANO_17_13, ROZ_SPD_USVIT_17_13, ROZ_CSSD_17_13, ROZ_KSCM_17_13), # názvy modelů
        file = "./ROZ_STRANY_2017_2013_HLM_OKRES.html", # output file
        custom.model.names = c("ANO 2017-13", "SPD 2017-13", "CSSD 2017-13", "KSČM 2017-13"), # název modelů
        reorder.coef = c( 2, 3, 4, 5, 6, 7, 8, 1),
        custom.coef.names = c("(Constant)","College degree (%)" , "Self-employed (%)", "Unemployed (%)" , "65+ years (%)", "Catholics (%)",
                              "Distraints (%)", "Number of inhabitants (log)"), # názvy proměnných
        custom.gof.names = c("AIC","BIC","Log Likelihood", 
                             "Num. obs.", "Num. groups: District",
                             "Var: District (Intercept)","Var: Residual"), 
        digits = 3, # decimal places for coefficients and SEs
        single.row = FALSE) # display BETA and SEs under each other

screenreg(list(ROZ_ANO_17_13, ROZ_SPD_USVIT_17_13, ROZ_CSSD_17_13, ROZ_KSCM_17_13),  # vizualizace modelů.... 
          digits = 3,
          custom.model.names = c("ANO 2017-13", "SPD 2017-13", "CSSD 2017-13", "KSČM 2017-13"))


performance::icc(NULL_ROZ_ANO_17_13)
performance::icc(NULL_ROZ_SPD_USVIT_17_13)
performance::icc(NULL_ROZ_CSSD_17_13)
performance::icc(NULL_KSCM_ROZDIL_17_13)





# Analýza ANO a SPD + doplnění komplet dat za KSČM a ČSSD - standardizované
sd_ROZ_ANO_17_13 <- lm(sd_ANO_ROZDIL_P ~ sd_Vzd_vysok_A + sd_Eko_podnikatel_A + sd_Eko_nezam_A + sd_Nepracduch_A + sd_Nab_rimsko_A + sd_EXEKUCE_OBCE + sd_INHABITANTS_2018_LOGGED, data=FINAL_17, weights=VOL_SEZNAM_17)
summary(sd_ROZ_ANO_17_13)

sd_ROZ_SPD_USVIT_17_13 <- lm(sd_SPD_USVIT_ROZDIL_P ~ sd_Vzd_vysok_A + sd_Eko_podnikatel_A + sd_Eko_nezam_A + sd_Nepracduch_A + sd_Nab_rimsko_A  + sd_EXEKUCE_OBCE + sd_INHABITANTS_2018_LOGGED, data=FINAL_17, weights=VOL_SEZNAM_17)
summary(sd_ROZ_SPD_USVIT_17_13)

sd_ROZ_CSSD_17_13 <- lm(sd_CSSD_ROZDIL_P ~ sd_Vzd_vysok_A + sd_Eko_podnikatel_A + sd_Eko_nezam_A + sd_Nepracduch_A + sd_Nab_rimsko_A + sd_EXEKUCE_OBCE + sd_INHABITANTS_2018_LOGGED, data=FINAL_17, weights=VOL_SEZNAM_17)
summary(sd_ROZ_CSSD_17_13)

sd_ROZ_KSCM_17_13 <- lm(sd_KSCM_ROZDIL_P ~ sd_Vzd_vysok_A + sd_Eko_podnikatel_A + sd_Eko_nezam_A + sd_Nepracduch_A + sd_Nab_rimsko_A + sd_EXEKUCE_OBCE + sd_INHABITANTS_2018_LOGGED, data=FINAL_17, weights=VOL_SEZNAM_17)
summary(sd_ROZ_KSCM_17_13)

library(texreg)
htmlreg(list(sd_ROZ_ANO_17_13, sd_ROZ_SPD_USVIT_17_13, sd_ROZ_CSSD_17_13, sd_ROZ_KSCM_17_13), # názvy modelů
        file = "./ROZ_STRANY_2017_2013_SD.html", # output file
        custom.model.names = c("ANO 2017-13", "SPD 2017-13", "CSSD 2017-13", "KSČM 2017-13"), # název modelů
        reorder.coef = c( 2, 3, 4, 5, 6, 7, 8, 1),
        custom.coef.names = c("(Konstanta)","Vysokoškoláci" , "Podnikatelé", "Nezaměstnaní" , "Důchodci", "Katolící",
                              "Distraints (%)", "Number of inhabitants (log)"), # názvy proměnných
        custom.gof.names = c("R2","Adj.R2","Počet pozorování", "RMSE"), # názvy proměnných
        digits = 3, # decimal places for coefficients and SEs
        single.row = FALSE) # display BETA and SEs under each other

#vytvoření grafu regresních koeficientů - standardizovaných
library(coefplot)
multiplot(sd_ROZ_ANO_17_13, sd_ROZ_SPD_USVIT_17_13, sd_ROZ_CSSD_17_13, sd_ROZ_KSCM_17_13, outerCI = 4, innerCI = 2, 
          intercept = FALSE,  title = "Změna volební podpory 2017 oproti 2013", xlab = "Hodnota koeficientu - standardizovaná  + 99% k.i.",
          ylab = "Název  koeficientu", sort = c("alphabetical"), pointSize = 3, color = "blue", cex = 5, 
          newNames=c(sd_Vzd_vysok_A="Vysokoškoláci",
                     sd_Eko_podnikatel_A="Podnikatelé",
                     sd_Eko_nezam_A="Nezaměstnaní",
                     sd_Nepracduch_A="Důchodci",
                     sd_Nab_rimsko_A="Katolící"),
          names=c(sd_ROZ_ANO_17_13= "ANO 2011", sd_ROZ_SPD_USVIT_17_13="SPD",
                  sd_ROZ_CSSD_17_13="ČSSD", sd_ROZ_KSCM_17_13="KSČM")) + scale_color_manual(values=c("blue", "orange", "red", "brown")) +
  geom_vline(xintercept=0, linetype="dashed", color = "black", size=1.1) + theme(plot.title = element_text(size=12,face="bold")) +
  theme(axis.text=element_text(size=13), axis.title=element_text(size=10,face="bold")) 

library(coefplot) #pouze podíl vysokoškoláků
multiplot(sd_ROZ_ANO_17_13, sd_ROZ_SPD_USVIT_17_13, sd_ROZ_CSSD_17_13, sd_ROZ_KSCM_17_13, outerCI = 4, innerCI = 2, 
          intercept = FALSE,  title = "Změna volební podpory 2017 oproti 2013", xlab = "Hodnota koeficientu - standardizovaná  + 99% k.i.",
          ylab = "Podíl vysokoškoláků v obci", sort = c("magnitude"), pointSize = 3, color = "blue", cex = 5, 
          newNames=c(sd_Vzd_vysok_A=" "), coefficients=("sd_Vzd_vysok_A"),
          names=c(sd_ROZ_ANO_17_13= "ANO 2011", sd_ROZ_SPD_USVIT_17_13="SPD",
                  sd_ROZ_CSSD_17_13="ČSSD", sd_ROZ_KSCM_17_13="KSČM")) + scale_color_manual(values=c("blue", "orange", "red", "brown"))



#porovnání analýzy sněmovních voleb 2017 s Prahou a bez Prahy
MODEL_CSSD_17 <- lm(PAR_2017_4_Českástr.sociálnědemokrat._p ~ Vzd_vysok_A + Eko_podnikatel_A + Eko_nezam_A + Nepracduch_A + Nab_rimsko_A, data=PS_KOMPLET_17_96, weights=PAR_2017_VOL_SEZNAM)
summary(MODEL_CSSD_17)

MODEL_CSSD_17_FULL <- lm(PAR_2017_4_CSSD_p ~ Vzd_vysok_A + Eko_podnikatel_A + Eko_nezam_A + Nepracduch_A + Nab_rimsko_A, data=FINAL_17, weights=VOL_SEZNAM_17)
summary(MODEL_CSSD_17_FULL)


screenreg(list(MODEL_CSSD_17, MODEL_CSSD_17_FULL),  # vizualizace modelů.... 
          digits = 3,
          custom.model.names = c("ČSSD 2017 BEZ PRAHY",
                                 "ČSSD 2017 S PRAHOU"))



# Hnutí ANO
ANO_17 <- lm(PAR_2017_21_ANO2011_p ~ Vzd_vysok_A + Eko_podnikatel_A + Eko_nezam_A + Nepracduch_A + Nab_rimsko_A, data=FINAL_17, weights=VOL_SEZNAM_17)
summary(ANO_17)

ANO_13 <- lm(PAR13_20_ANO_2011_p ~ Vzd_vysok_A + Eko_podnikatel_A + Eko_nezam_A + Nepracduch_A + Nab_rimsko_A, data=FINAL_17, weights=VOL_SEZNAM_17)
summary(ANO_13)



screenreg(list(ANO_17, ANO_13),  # vizualizace modelů.... 
          digits = 3,
          custom.model.names = c("ANO 2011 PS2017",
                                 "ANO 2011 PS2013"))

library(texreg)
htmlreg(list(ANO_17, ANO_13), # názvy modelů
        file = "./MODELY_ANO.html", # output file
        custom.model.names = c("ANO 2011 PS2017",  "ANO 2011 PS2013"), # název modelů
        reorder.coef = c( 2, 3, 4, 5, 6,1),
        custom.coef.names = c("(Konstanta)","Vysokoškoláci" , "Podnikatelé", "Nezaměstnaní" , "Důchodci", "Katolící"), # názvy proměnných
        custom.gof.names = c("R2","Adj.R2","Počet pozorování", "RMSE"), # názvy proměnných
        digits = 3, # decimal places for coefficients and SEs
        single.row = FALSE) # display BETA and SEs under each other



# Model Tomáš
#####################################################################################################

FINAL_17$POPULISTE_ROZDIL <- (FINAL_17$ANO_ROZDIL_P + FINAL_17$SPD_USVIT_ROZDIL_P)


MODEL1 <- lm(POPULISTE_ROZDIL ~ CSSD_ROZDIL_P + KSCM_ROZDIL_P + ODS_ROZDIL_P, data=FINAL_17)
summary(MODEL1)


MODEL2 <- lm(POPULISTE_ROZDIL ~ Vzd_vysok_A + Eko_podnikatel_A + Eko_nezam_A + Nepracduch_A + Nab_rimsko_A, data=FINAL_17)
summary(MODEL2)

MODEL3 <- lm(POPULISTE_ROZDIL ~ CSSD_ROZDIL_P + KSCM_ROZDIL_P + ODS_ROZDIL_P + Vzd_vysok_A + Eko_podnikatel_A + Eko_nezam_A + Nepracduch_A + Nab_rimsko_A, data=FINAL_17)
summary(MODEL3)

library(texreg)
htmlreg(list(MODEL1, MODEL2, MODEL3), # názvy modelů
        file = "./MODELY_ROZDILY.html", # output file
        custom.model.names = c("Model 1",  "Model 2", "Model 3"), # název modelů
        reorder.coef = c( 2, 3, 4, 5, 6, 7, 8, 9, 1),
        custom.coef.names = c("(Konstanta)","ČSSD rozdíl (%)", "KSČM rozdíl (%)", "ODS rozdíl (%)","Vysokoškoláci" , "Podnikatelé", "Nezaměstnaní" , "Důchodci", "Katolící"), # názvy proměnných
        custom.gof.names = c("R2","Adj.R2","Počet pozorování", "RMSE"), 
        digits = 3, # decimal places for coefficients and SEs
        single.row = FALSE) # display BETA and SEs under each other


MODEL3 <- lm(POPULISTE_ROZDIL ~ CSSD_ROZDIL_P + KSCM_ROZDIL_P + ODS_ROZDIL_P + Vzd_vysok_A + Eko_podnikatel_A + Eko_nezam_A + Nepracduch_A + Nab_rimsko_A  + EXEKUCE_OBCE + INHABITANTS_2018_LOGGED, data=FINAL_17)
summary(MODEL3)


#### v ČLÁNKU!!! table 1!
library(texreg)
htmlreg(list(ROZ_ANO_17_13, ROZ_SPD_USVIT_17_13, ROZ_CSSD_17_13, ROZ_KSCM_17_13, MODEL3), # názvy modelů
        file = "./TABLE1.html", # output file
        custom.model.names = c("ANO 2017-13", "SPD 2017-13", "CSSD 2017-13", "KSČM 2017-13", "ANO+SPD 17-13"), # název modelů
        reorder.coef = c( 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 1),
        custom.coef.names = c("(Constant)","College degree (%)" , "Self-employed (%)", "Unemployed (%)" , "65+ years (%)", "Catholics (%)",
                              "Distraints (%)", "Number of inhabitants (log)",
                              "ČSSD change", "KSČM change", "ODS change"), # názvy proměnných
        custom.gof.names = c("R2","Adj.R2","Number of cases", "RMSE"), # názvy proměnných
        digits = 3, # decimal places for coefficients and SEs
        single.row = FALSE) # display BETA and SEs under each other








# S vážením
MODEL1 <- lm(POPULISTE_ROZDIL ~ CSSD_ROZDIL_P + KSCM_ROZDIL_P + ODS_ROZDIL_P, data=FINAL_17, weights=VOL_SEZNAM_17)
summary(MODEL1)

MODEL2 <- lm(POPULISTE_ROZDIL ~ Vzd_vysok_A + Eko_podnikatel_A + Eko_nezam_A + Nepracduch_A + Nab_rimsko_A, data=FINAL_17, weights=VOL_SEZNAM_17)
summary(MODEL2)

MODEL3 <- lm(POPULISTE_ROZDIL ~ CSSD_ROZDIL_P + KSCM_ROZDIL_P + ODS_ROZDIL_P + Vzd_vysok_A + Eko_podnikatel_A + Eko_nezam_A + Nepracduch_A + Nab_rimsko_A, data=FINAL_17, weights=VOL_SEZNAM_17)
summary(MODEL3)



#####################################################################################################
################################# Korelační analýza #################################################
#####################################################################################################



#KORELACNI MATICE 2017
PS_corr17 <- PS_KOMPLET_17_96[, c("PAR_2017_21_ANO2011_p", "PAR_2017_1_Občanskádemokratickástrana_p", "PAR_2017_15_Českápirátskástrana_p", "PAR_2017_8_Komunistickástr.ČechaMoravy_p", "PAR_2017_4_Českástr.sociálnědemokrat._p",
                                   "PAR_2017_29_Svob.apř.dem.T.OkamuraSPD_p", "PAR_2017_20_TOP09_p", "PAR_2017_24_Křesť.demokr.unieČs.str.lid._p", "PAR_2017_7_STAROSTOVÉANEZÁVISLÍ_p",
                                   "Nab_rimsko_A", "Vzd_vysok_A", "Eko_podnikatel_A", "Eko_nezam_A", "Nepracduch_A")]

PS_korelace_17 <-cor(PS_corr17, use="complete.obs", method="pearson")
colnames(PS_korelace_17) <- c("ANO PS 2017", "ODS 17", "Piráti 17", "KSČM 17",
                               "ČSSD 17","SPD 17", "TOP09 17", "KDU-ČSL 17", "STAN 17", "Katolící", "Vysokoškoláci" , "Podnikatelé", "Nezaměstnaní" , "Důchodci")

rownames(PS_korelace_17) <- c("ANO PS 2017", "ODS 17", "Piráti 17", "KSČM 17",
                               "ČSSD 17","SPD 17", "TOP09 17", "KDU-ČSL 17", "STAN 17", "Katolící", "Vysokoškoláci" , "Podnikatelé", "Nezaměstnaní" , "Důchodci")

library(corrplot)
corrplot(PS_korelace_17, method="number", type="upper", tl.col = "black",  number.cex = 0.8, cl.cex = 0.8, tl.cex = 0.7)



#KORELACNI MATICE ANO, ČSSD, KSČM, SPD
#######################################
PS_corr_LEVICE <- PS_KOMPLET_17_96[, c("PAR_2017_21_ANO2011_p", "PAR_2017_4_Českástr.sociálnědemokrat._p", "PAR_2017_8_Komunistickástr.ČechaMoravy_p", 
                                  "PAR_2017_29_Svob.apř.dem.T.OkamuraSPD_p",
                                  "PAR13_20_ANO_2011_p", "PAR13_1_Česká_strsociálně_demokrat_p", "PAR13_21_Komunistická_strČech_a_Moravy_p",
                                  "PAR13_17_Úsvit_přímé_demokrTOkamury_p",
                                  "PAR10_9_Česká_strsociálně_demokrat_p", "PAR10_6_Komunistická_strČech_a_Moravy_p", "PAR10_4_Věci_veřejné_p",
                                  "PAR06_10_Česká_strsociálně_demokrat_p", "PAR06_20_Komunistická_strČech_a_Moravy_p",
                                  "PAR02_3_Česká_strsociálně_demokrat_p", "PAR02_23_Komunistická_strČech_a_Moravy_p", 
                                  "PAR98_7_Česká_strsociálně_demokrat_p", "PAR98_9_Komunistická_strČech_a_Moravy_p", 
                                  "PAR96_2_Česká_strsociáldemokratická_p", "PAR96_15_Komunistická_strČech_a_Moravy_p", 
                                  "PAR96_17_Sdružpro_repRepublstrČsl_p", 
                                  "Nab_rimsko_A", "Vzd_vysok_A", "Eko_podnikatel_A", "Eko_nezam_A", "Nepracduch_A")]

PS_korelace_LEVICE <-cor(PS_corr_LEVICE, use="complete.obs", method="pearson")
colnames(PS_korelace_LEVICE) <- c("ANO 17","ČSSD 17", "KSČM 17", "SPD 17",
                                  "ANO 13","ČSSD 13", "KSČM 13", "ÚSVIT 13",
                                  "ČSSD 10", "KSČM 10", "VV 10",
                                  "ČSSD 06", "KSČM 06", "ČSSD 02", "KSČM 02", 
                                  "ČSSD 98", "KSČM 98", "ČSSD 96", "KSČM 96", "SPR-RSČ 96", 
                                  "Katolící", "Vysokoškoláci" , "Podnikatelé", "Nezaměstnaní" , "Důchodci")

rownames(PS_korelace_LEVICE) <- c("ANO 17","ČSSD 17", "KSČM 17", "SPD 17",
                                  "ANO 13","ČSSD 13", "KSČM 13", "ÚSVIT 13",
                                  "ČSSD 10", "KSČM 10", "VV 10",
                                  "ČSSD 06", "KSČM 06", "ČSSD 02", "KSČM 02", 
                                  "ČSSD 98", "KSČM 98", "ČSSD 96", "KSČM 96", "SPR-RSČ 96", 
                                  "Katolící", "Vysokoškoláci" , "Podnikatelé", "Nezaměstnaní" , "Důchodci")
library(corrplot)
corrplot(PS_korelace_LEVICE, method="number", type="upper", tl.col = "black",  number.cex = 0.8, cl.cex = 0.8, tl.cex = 0.7)


# Korelogramy výsledků
object <- ggplot(PS_KOMPLET_17_96, aes(x = PAR10_9_Česká_strsociálně_demokrat_p, y = PAR_2017_21_ANO2011_p, color = Kraj_nazev)) +
  geom_point() + geom_smooth(method=lm, fullrange=TRUE, se=FALSE, color = "black") + labs(x = "ČSSD podíl hlasů 2010", y = "Zisk ANO 2017")
library(plotly)
ggplotly(object)


plot_ly(PS_KOMPLET_17_96, x = ~PAR10_9_Česká_strsociálně_demokrat_p, y = ~PAR_2017_21_ANO2011_p, type = 'scatter', mode = 'markers',
        color = ~Kraj_nazev,hoverinfo = 'text',
        text = ~paste('</br> Kraj: ', Kraj_nazev,
                      '</br> Název obce: ', Nazev_obce,
                      '</br> Okrsek: ', OKRSEK_16,
                      '</br> ANO 2011: ', PAR_2017_21_ANO2011_p,
                      '</br> ČSSD 2010: ', PAR10_9_Česká_strsociálně_demokrat_p)) %>% 
  layout(title = 'Korelogram ANO 17 a ČSSD 10', xaxis = list(title = 'Zisk ČSSD ve sněmovních volbách 2010'),
         yaxis = list(title = 'Zisk ANO ve sněmovních volbách 2017'))




plot_ly(PS_KOMPLET_17_96, x = ~PAR13_1_Česká_strsociálně_demokrat_p, y = ~PAR_2017_4_Českástr.sociálnědemokrat._p, type = 'scatter', mode = 'markers',
        color = ~Kraj_nazev,
        text = ~paste('Název obce: ', Nazev_obce, '- okrsek: ', OKRSEK_16)) %>% layout(title = 'Korelogram ČSSD', xaxis = list(title = 'Zisk ČSSD ve sněmovních volbách 2013'),
                                                                                       yaxis = list(title = 'Zisk ČSSD ve sněmovních volbách 2017'))

plot_ly(PS_KOMPLET_17_96, x = ~PAR13_1_Česká_strsociálně_demokrat_p, y = ~PAR_2017_4_Českástr.sociálnědemokrat._p, type = 'scatter', mode = 'markers',
        color = ~Kraj_nazev, hoverinfo = 'text',
        text = ~paste('</br> Název obce: ', Nazev_obce,
                      '</br> Okrsek: ', OKRSEK_16,
                      '</br> ČSSD 2017: ', PAR_2017_4_Českástr.sociálnědemokrat._p,
                      '</br> ČSSD 2013: ', PAR13_1_Česká_strsociálně_demokrat_p)) %>% 
  layout(title = 'Korelogram ČSSD', xaxis = list(title = 'Zisk ČSSD ve sněmovních volbách 2013'),
         yaxis = list(title = 'Zisk ČSSD ve sněmovních volbách 2017'))

